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Abstract 



In the present work the cohapse scenario of some exact non-spherical models with a 
CJ ' minimally coupled scalar field is studied. Scalar field collapse with planar as well as toroidal, 

cylindrical and pseudoplanar symmetries have been investigated. It is shown that the scalar 
field may have collapsing modes even if it has the equation of state corresponding to that of 
a dark energy. 



> 

in 1 Introduction 

\o ■ 

•r^ ', It is now well known that a collapsing matter distribution may not always lead to the formation 
i/-^ [ of a black hole. The spacetime singularity, formed out of the unhindered collapse, may not 
O ' be covered by an event horizon and thus exposed to an external observer. The formation of 
^ ■ black holes, for a spherically symmetric dust distribution, normally depends on the initial density 
, proffie[l]. A star is expected to collapse under its own gravity when its nuclear fuel gets exhausted 
so that there is no pressure to balance gravity. So the approximation of the matter distribution 
as dust is quite a well-motivated one. But as there is not much phenomenological example of 
the equation of state of the high density collapsing fluid, investigations regarding the collapse of 
other forms of matter become important. Scalar fields naturally step in as a candidate for other 
forms of matter, mainly because of their ability to mimic different kinds of matter depending on 
the choice of potentials. A lot of spherical scalar field models already figure in the literature. 
Some of them involve exact analytical study [2] and some are numerical endeavours [3]. Careful 
investigations regarding a scalar field collapse bring out the fascinating possibility of a critical 
phenomenon [4] . 

Although mathematically elegant, a spherically symmetric collapse, is hardly found in nature 
as most of the stellar bodies do show departure from sphericity. One may approach the problem 
in two different ways. One is to look at exact non-spherical models and the other resorts to 
a perturbative approach where a departure from the spherical symmetry in collapse scenario is 
considered. For an extensive review, we refer to the work by Miller and Sciama[5]. 

The present work investigates some exact non-spherical models with a minimally coupled scalar 
field. Bronnikov and Kovalchuk made an extensive study of the various non-spherical collapse 
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with a pressureless dust distribution [6] and then generahzed the same to include a charged dust in 
planar symmetry [7] and later an electromagnetic field for various non-spherical distributions [8]. 
The present paper extends the work for scalar fields in planar, cylindrical, toroidal and pseudo- 
planar symmetries. 

As the interest in scalar fields has increased manyfold because of their role as a possible 
candidate as a 'dark energy' driving the present accelerated expansion of the Universe(see ref[9] 
for reviews), the study of the collapse of a scalar field distribution becomes important as this 
may lead to some idea about the clustering properties of the dark energy which are yet to be 
ascertained. 

In section 2 we write down the Einstein's equations for the scalar field distribution for a general 
kind of metric which gives rise to different symmetries as special cases. Section 3 discusses the 
collapsing models for a planar symmetry, section 4 describes scalar field collapse for cylindrical, 
toroidal and pseudoplanar symmetry. In section 5 we make some concluding remarks. 

2 Field Equations 

We take a metric. 



where 7, A, a and ui depend on R and T. For different choices oi f — f{r)) and ou, one has different 
symmetries as, 

/ = sin77, uj — (spherical), 
/ = sinh?7, UJ — (pseudospherical) , 
f — 1, oj ^ (pseudoplanar, cylindrical, toroidal), 
/ = 1, a; = (planar). 

In the present work the collapse scenario for the latter two are investigated. The energy mo- 
mentum tensor T^^' for scalar field where is a function of 'i?' and T'and T^(0) is the scalar 
potential is given as, 

]g"^cj>,acl>,p - V^(0)1 g^.u (1) 



The Einstein equations for the said metric with T^^, given as in (1) with f — 1 are 



27-e-^^ - e-^^(2^ + %+^') + 2y-e-^^ + (^ - J^)e-^>^ = J(e-^V' + e-^^0^) - ^(0), (2) 
a a a a'^ 2 



-2t(-A - - f - 2f a; - d;2 - a) + 7A + 7f + 70; - Af - Aa;)+ 

-2A(^+2-a;'W'W'+7''+7''-A'--AV-7'A'+7'-+7V) = ^(-e-2V'+e-V)-V^(0), 
a a a a 2 



e-2T(-A - A^ - f + 2^0; - a;2 + c;} + 7A + 7! - 7^^ - A^ + Xuj)+ 



(3) 



(4) 



a a a a 2 
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e-'\2^ + ^ + o;'^ - 2A'-) - e-'^{2\^ + ^ " = ^(-e-'V" - e'^^^) - V^(0), (5) 
a a a a'^ I 

-2(-+a;'a;-A--7-) = </>0'. (6) 
a a a 

where we have taken SnG = 1, a dot and a prime represent differentiations with respect to time 
'T' and space coordinate 'R' respectively. Equations (2) — (6) represent Gi^, G^, Gq^ and 
Gi^ equations repectively. 

In order to make the equation system tractable, we assume ^{R,T), \{R,T) and a{R,T) 
are all separable as products of functions of R and that of T. The T-dependence of 7 and in- 
dependence of A can now be absorbed in dT"^ and dR^ respectively by suitable scaling of T and 
R without any further loss of generality. Now one has, 

7 = 7(i?), 

A = A(T), 
a{R,T) = a{R)l3{T). 



3 Planar symmetry (cj = 0, / = 1) 

For u! — the field equations take the form. 



e-"(< + ^)--^4.})4e-V-VW. (7) 

+ y + + j'^) + e-'^(-X - - I - A^) = le-'-'4,' - (8) 
a a p p 2 



a p p 



(2— + ^) - e-^'(2A^ + ^) = --e-^V - V(0), (9) 



4^ + 2A^ + 2y| = 0. (10) 

pa a p 



Integrating equation(lO) one gets, 

A 7 



In /9 In a ' 

where a constant of integration is put equal to zero. From equations (8) and (9) we get. 



VW ^ + ^ + 7" + + - i.-(-A - A= - 3A^ - I - I). (12) 



and 



(11) 



4>' = e-\^-(-^ - 4 + 7" + + y-) - A - A^ - I + a| + ^. (13) 
q; a p p p^ 

In the above system of equations is assumed to be a function of T alone. This is a simplifying 
assumption made to solve the system of equations easily. Therefore what follows from the above 
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assumption is that 1^(0) is also a function of T and to implement this, we put 7 = constant. The 
value of the constant is taken as zero without any further loss of generality. Equations (12) and 
(13) then give, 

a a , , 
3— + ^ = A;i, 14 

II /2 

^ = -h, 15 

where ki and k2 are constants. Both the equations (14) and (15)are satisfied by the equation 

a = e^^, (16) 

and ki and ^2 are related as ki = 2k2 for the sake of consistency. The above form of the solution 
of a{R) also satisfies the other equations (7) - (10). With 7 = 0, equation (11) therefore has the 
form 

A = ln^. (17) 
With this the metric takes the form ds^ ^ ^^^7^2 _ ^2^^2 _ ^2g^ii^^^2 ^ ^^2^ g^j^ 
find 13. 

We make an assumption that the potential is directly proportional to the kinetic energy term, 

i.e., 

V{(l)) = m(/)^ (18) 

where m is a constant. As goo = 1 and = 0(t), 0^ is actually which is indeed a scalar 

and equation (18) is thus consistent. However, it deserves mention that the field in this case is 
not a k-essence field where the Lagrangian is given by L = L{X) where X — 0^^^''^. With the 
choice of 7 = 0, the energy momentum tensor for present model looks like 

Ti; = {(j)fdiag {{1/2 + m), (1/2 - m), (1/2 - m), (1/2 - m)). 

Actually the model does have a potential, and the present work deals with those which inci- 
dentally are proportional to the kinetic part. For instance, the simple case where is spatially 
homogeneous, a choice like 

y(0) oc exp{-^) 

will in fact be a possibility that gives rise to equation (18) at least in a special case. 

The equation of state parameter, w = ^ = i|^2_|_y(^) "^i^ be a constant given hy w = 
For any finite positive value of m, w remains greater than —1 and hence the scalar field is in fact 
a quintessence field. 

Putting the expressions of T^(0) and 0^ from (12) and (13) one gets, 

P'-^^ + Csr'^\ (19) 

where Ci = , C2 = 2(i+2mj ^^"^ ^3 integration constant. Now since it becomes difficult 

to further integrate equation (19) to obtain /3{T), we solve it numerically. From (19) we have. 



-y'^ + Car^^^ (20) 
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As we are interested in a collapsing model, the negative root is considered. f3{T) is plotted as 
a function of T for different values of m. For m = i.e. a scalar field model without any 
potential, we have /3 = —\J^ + Csf3~'^ and giving some appropriate values of the constants ki 
and C3 and applying reasonable boundary conditions it is seen from (figure 1) that (3{T) starts 
from a finite value but becomes zero at a finite time. As the proper volume is proportional to 
{y/—g = this indicates that the distribution collapses to the singularity of a zero proper 

volume. Different reasonable sets of values for ki and C3 would lead to different values for time 
T when the singularity is attained but the nature of the collapse remains similar. 



For m = 




0.0 0.1 0.2 0.3 0.4 

>T 



Figure 1: Plot of /3(T) vs T for m = with ki = 3 and C3 = 2 (w = 0) 
for different values of the constant (figure 2 shows one of them). 




Figure 2: Plot of /3{T) vs T for m = 1/4 with ki = 1 and C3 = 1 {u = 0) 

For these collapsing model is not surprising as the matter content satisfies the energy 

conditions (p + 3p > 0). If m > 1/2, the effective pressure becomes negative and when m > 1, 
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the energy condition is violated and the scalar field is apt to act as a dark energy. As an example 
we choose m = 3/2, for which equation (19) yields /3 = —\J^ + C^P^/"^. With some reasonable 
boundary conditions and putting the values of constants ki and C3 both equal to 1 one can get 
a collapsing model as seen in (figure 3). 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

>T 

Figure 3: Plot of /3(T) vs T for m = 3/2 with ki = 1 and C3 = 1 (a; = 0) 

For values of m less than 1/2 we have sharply collapsing models where the volume falls to 
zero very rapidly particularly towards the end. For m = 3/2, the rate of collapse slightly fiattens 
out towards the end, but one can have a collapsing model none-the-less. 

It is relevant to point out that at the singularity (/3(T) — )■ 0), both 0^ and Vi^cj)) would blow 
up leading to an infinite energy density and pressure. Equation (13) now yields 

= (0ir' + 02r")^ 

where 02 are positive constants, related with the constants of integration and m, and n = 
6/(1 + 2m). This cannot be analytically integrated for all values of m. If one chooses m = 1, it 
can be shown that oc InT. Using the ansatz (18), one can now easily show that 

V{(f)) = me"^"^. 

4 Cylindrical, toroidal, pseudoplanar symmetry {uu ^ 0) 

Now some models with a non-zero u is considered which would lead to cylindrical, toroidal or 
pseudoplanar kind of symmetry depending on the range of the coordinates. When the coordinates 
T) and grange from — oo to +00, one has a pseudoplanar symmetry. If one of the coordinates rj or 
^ is an angle coordinate ranging from to 27r, one would have a cylindrical symmetry. If both of 
them are angles, one has a toroidal symmetry. In what follows, we assume u is either a function 
of R or that of T for the sake of simplicity. 
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4.1 u = uj{R) 

With CO = U!{R) and the same simpUfications as before, that is, 7 = j{R), A — A(T), a{R,T) — 
a{R)/3{T) and = 0(7") the field equations are as follows, 

+ ^ - - + {) - ie-^^ - V(,). (21) 

e-"(^+2u;'-+u;''W+7"+7"+7 -+T'i^')+e"'"(-A-A"-f -A^) = I- e-'^i'--V {,!,), (22) 
a a ct p p 2 

q; q; a p p 2 

+ ^ + '^'^) - e-(2Af + {) ^ -ie->^ - n^), (24) 



ck' _ , . a' 



2uj'uj + 2X— + 27'^ = 0. (25) 
a a p 



Since uj = U!{R) only, equation (25) takes the form. 



-2^- + 2A- + 27'^ = 0. (26) 
pa a p 



Integrating this equation we get 

A 7 



ln/3 ^ Ina' ^^^^ 

exactly the same relation as in equation (11). By the same logic of being a function of T only, 
7 is a constant and is chosen to be zero. Equations (22) and (23) yield an equation of the form 

^ + 2^ = 0, (28) 

which on integrating gives, 

^' = 4' (29) 

where ^ is a constant of integration. Substituiting ou' from (29) in the field equations, one has an 
equation of the form, 

^-- = K, (30) 
where i^T is a constant. Integrating the above equation one gets the solution of a{R) as 



a(R) = Ai^sinh(i? + C) (31) 

where Ai = and C is an integration constant and the solution is consistent with all the field 
equations. Now putting the solution of a{R) from (31) and integrating equation (29) one obtains 

u^Kln (te^nh • (32) 
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Using equations (21) and (24) we have 



^^^2^-2^-2^. 

Using equations (33), (34) and the ansatz (18) one can write, 

•2 K{l-2m) ^^^1(1^ 

P = -7^. r + ^" 

2(1 — m) 



(33) 
(34) 



(35) 



where E is an integration constant. Taking the positive root of $ gives an ever expanding modeL 
But considering the negative root there is a possibihty to have a collapsing model and hence the 
negative root is chosen. 

Putting the value of m = and constants E and K equal to 1, equation (35) looks like 
/3 = —'fj Y + EI3~'^. From the numerical plot of f3{T) vs T in (figure 4) it is seen that the volume 
starts from non-zero value at some initial time T = and then becomes zero at a finite value 
of time T. Putting some other values of the constants only shifts the time at which the proper 
volume crushes to the singularity but the nature of the graphs remain unchanged. 



A value of m = 1/4 that ie V{(j)) 



1^2 



and same values of the constants also yield a model 



A 




->T 



Figure 4: Plot of f3{T) vs T for m = with K = 1 and E = 1 {u = u{R)) 

collapsing in nature as is evident from (figure 5). 
As before we take same example of m = 3/2 to represent a dark energy candidate. Equation 

(35) now becomes /3 = —\J2K + E/d^/"^. It is evident from the plot in (figure 6) that the volume 
becomes zero gradually at some value of T. In this case also the different values of the constants 
only change the time for which (3{T) keeping the nature of the graphs unchanged. 



Equation (34) and (35) yield 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 

>T 

Figure 5: Plot of /3(T) vs T for m = 1/4 with i^T = 1 and ^ = 1 (w = uj{R)) 




Figure 6: Plot of (3{T) vs T for m = 3/2 with K = A and E = 2 (oj = uj{R)) 

= (03r' + 04r')^ 

where 03, 04 are constants, related with the constants of integration and m, and / = 6/(1 + 2m). 
This cannot in general be integrated. But (p indeed diverges when /3, and hence the proper volume, 
goes to zero. Thus the energy denisty has a singularity when /3 goes to zero. 

4.2 CO = u{T) 

Taking u = uj{T) the field equations can be written as, 

e-'H^y- + 4) - ^-"'(4 + ^ + = l'"""^" - ^^'^)' ^36) 
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e-^^(- + 7" + yVy-) + e-^-(-A-A^-^-2^a;-a;^-d;-A^-Aa;) = le'V-V^l^), (37) 
a a p p pi 

e-'\— + 7" + y ' + 7 -)e-'^(-A - - f + 2^(1; - (0^ + a; - a| + A(:;) = ^e'^^./.^ - 1^(0) , (38) 
a a p p p 2 

-2|^ + 2A- + 27'| = 0, (40) 
pa a p 

the same simplifying assumptions like 7 = 7(-R), A = A(T), a{R,T) = a{R)/3{T) and = 0(T) 
are kept intact. Following similar arguments as before, one can write, A = In ^. Solving the above 
field equations (36) to (39) one easily arrives at 

a = e^^, (41) 

which is consistent with all the field equations, Ki being a constant of integration. Now subtract- 
ing equation (38) from (37) one obtains after integration, 

- = ^' (42) 
where S is a constant of integration. Equations (36) and (39) give, 

m = (§ + 2^)-e-^(^ + ^), (43) 
p P'^ a 



and 



^--2,^H-.^H-^)-2e-< (44) 



Use of equations (17), (41) and (42) in equations (43) and (44) results in 



-A' = -2(^ + « + ^) - 27^- («) 



and 

J 52 /32 
- H V —)- 2— 

Making the same assumption as in equation (18) and integrating a combination of equations (45) 
and (46) one gets, 

(3^ = _ B ^^-4(l+n^) (47) 

^ 2(l + m) 2^ ^ ' ^ ' 

where D is a constant of integration. Choosing the negative root of /9(T), the plots of /9(T) vs T 
reveal that for m = and m— 1/4, with appropriate values of the constants, the volume crushes 
to zero sharply at some finite time (figure 7 and figure 8). Plotting I3{T) vs T for m = 3/2 
also, with suitable value of the constants, one gets a collapsing model (figure 9). The point to be 
noted is that in the previous models for a; = or a; = uj{R), m = 3/2 i.e., when the scalar field 
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0.0 ^ i\ 

0.000 0.005 0.010 0.015 0.020 0.025 
>T 

Figure 7: Plot of (3{T) vs T for m = with Ki = 0.002, B = 9 and D = 150 (w = w(T)) 




0.2 - 



0.0 b , , , , \ , , , , \ , , , , \ , , , , \ , , , , \ , , , , \ , , , J 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 

>T 

Figure 8: Plot of /3(T) vs T for m = 1/4 with Ki = Q, B = 9 and D = 10 {oo = uj{T)) 

is consistent with a dark energy candidate, we had collapsing models where the volume becomes 
zero in a sedate manner, but in this case when u = uj{T) even for values of m > 1/2 we have 
collapsing models where the volume falls to zero quite abruptly. 

In this case also one can write as a pwer function of /3 as 

where 05, 06 and 07 are constants, related with the constants of integration and m. One cannot 
intergrate for except for very special cases. Like the other cases, one still has a divergence of 
energy density at the singularity of zero proper volume. 
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1.0 
0.9 

j 0.7 
I 0.6 
0.5 



0.000 0.005 0.010 0.015 
>T 

Figure 9: Plot of (3{T) vs T for m = 3/2 with Ki = -0.002, B = 0.2 and D = 100 (w = uj{T)) 

5 Discussions 

Non spherically symmetric matter distributions with a minimally coupled scalar field can lead to 
collapsing models. With some choices of the arbitrary constants of integration, the models evolve 
to produce singularity of a zero proper volume at some finite future. The intriguing feature is 
that even for the cases where the scalar field acts as a dark energy candidate, one can have 
collapsing solutions with the only difference that for these cases the volume gradually becomes 
zero unlike the cases where m < 1/2 which leads to a suddenly collapsing situation at a finite 
future. However, in one case, u = uj{T), even the dark energy candidates may lead to a sudden 
collapse. For a positive value of m, it is obvious that w > —1, that is, the model will never be a 
phantom. It infact mimics a quintessence field. 

This work is indeed a specal case where the potential is assumed to be proportional to the 
kinetic part, but actually the model does not mimic the k-essence model. Here one does have a 
potential. At least in one special case, for the planar symmetry, an example of the potential is 
given which can give rise to this kind of a behaviour. 

This work also clearly indicates that the clustering properties of a 'dark energy' indeed deserves 
a more critical attention. 
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